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TECHNICAL NOTE 2367 


GENERAL PLASTIC BEHAVIOR AND APPROXIMATE SOLUTIONS OF 
ROTATING DISK IN STRAIN-HARDENING RANGE 
By M. H. Lee Wu 


SUMMARY 

A partly linearized solution of plastic deformation of a 
rotating disk “based on the deformation theory of plasticity and 
considering finite strains is obtained. The stresses and the 
strains of this problem for a given material and a given maximum 
strain can be obtained merely by a simple multiplication using the 
tables presented herein. This method is used to investigate the 
general plastic behavior of a rotating disk. An approximate method 
is also given in which the stresses are calculated by using strains 
obtained from the ideally plastic material, and the tensile true 
stress -strain curve of the material. 

Numerical examples are calculated by the two methods and 
agree very well with the exact solution based on deformation the- 
ory previously obtained. Calculations are also made for ideally 
plastic material and for the power-function approximation for pur- 
poses of comparison. 

The following conclusions, similar to those resulting from 
the linearized solution of the thin plate with a circular hole, are 
obtained for this problem: 

(1) The variation of a parameter, which is determined from the 
octahedral shear stress-strain curve of the material, can be used 
as a general criterion of the applicability of deformation theory. 

(2) The ratios of strain along the radius to the maximum value 
and the ratios of principal stresses are essentially independent 

of the octahedral shear stress-strain curve of the material, but 
the distributions of the stresses, and therefore the rotating 
speeds of the disk, depend very much on the material. 

(3) The results obtained from the ideally plastic material 
with the infinitesimal strain concept give good approximate values 
of strains but not of stresses. 
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(4) The rotating speed of the disk for a given maximum strain 
of the disk can he determined directly from the tensile true stress - 
strain curve of the material. 

(5) Good correlation between the experimentally determined 
bursting speed and the calculated value determined directly from 
the tensile stress-strain curve of the material is obtained. 

(6) If a simple analytical function representing the octa- 
hedral shear stress -strain relation is required for analysis, the 
power-law approximation can be used. 


INTRODUCTION 

In the design of a high-speed or a highly stressed machine 
member such as a turbine or a compressor rotor, the distributions 
of stresses and strains in the strain-hardening range must be 
known. The problem of a rotating disk for an ideally plastic 
material was solved by Nadai (reference l) . A solution for gas- 
turbine disks considering small plastic strain in the strain- 
hardening range is obtained in reference 2 by a trial-and-error 
procedure and by using elastic stress and strain distributions as 
a first approximation. An experimental investigation for the high- 
speed rotating disk is made in reference 3; distributions of plas- 
tic strains (logarithmic strains) for different types of disk are 
measured. The effect of strength and ductility on the burst char- 
acteristics of rotating disks are experimentally investigated in 
reference 4. An exact solution based on deformation theory for 
plane plastic stress problems with axial symmetry (including a cir- 
cular membrane under pressure, a rotating disk, and a thin plate 
with a circular hole) in the strain-hardening range is obtained 
in reference 5j numerical calculations are made for Inconel X and 
Timken alloy 16-25-6. A linearized solution of plastic deformation 
of a thin plate with a circular hole is given in reference 6. This 
linearized solution is not only simple and accurate but also offers 
a means of investigating the general plastic behavior of that prob- 
lem for most materials . 

Extension of this method to the problem having the additional 
complication of body forces, such as a rotating disk, is therefore 
interesting. (For a circular membrane under pressure, the normal 
pressure can be treated in a manner similar to the centrifugal 
force of a rotating disk.) The partly linearized solution obtained 
at the NACA Lewis laboratory and presented herein is also used to 
investigate general plastic behavior for this problem. 
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SYMBOLS 

The following symbols are used in this report: 


A^3,C, 
> D,E,F 

coefficients of nonlinear differential equations; functions 
of a, x } and r/k 

w 

b 

original outer radius of rotating disk 

H,J,L 

trigonometric functions of a 

h 

instantaneous thickness of disk 

h. .. 
mit 

initial thickness of disk 

K 

strain-hardening constant 

K 1 

arbitrary loading constant 

k 

constant, in dimension of length 

m 

parameter relating to strain hardening 

n 

parameter relating to criterion of applicability of deforma- 
tion theory 

r 

radial coordinate of undeformed disk 

u 

radial displacement 

z 

axial coordinate 

a 

parameter indicating ratio of principal stresses 

r 

octahedral shear strain 

€ 

logarithmic strain (natural strain), logarithm of instantaneous 
length divided by initial length of element 

e 

angular coordinate 

p 

mass per unit volume 

a 

normal true stress, force per unit instantaneous area 

T 

octahedral shear stress 

- CD 

angular velocity 
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Subscripts : 

b at outer radius b 

o at center of disk 

r,0,z principal directions; radial, tangential, and axial directions 


BASIC EQUATIONS 

A disk having an original outer radius b and an initial thick- 
ness h^ n i-k rotating about its axis with an angular velocity cd 
is shown in figure 1(a) . A small element defined by A0 and 
A(r+u) taken at the radius (r+u) in the deformed state is given 
in figure l(b) . In the undeformed state, this element is located 
at r and defined by A 9 and Ar. The instantaneous thickness h 
of the element and the stresses acting on the element are also 
shown in the figure. 

The relations of stresses and strains based on the deformation 
theory for plane plastic stress in the cylindrical coordinates are 
(references 7 and 8): 


e r +£ e + € z = 0 

(1) 

T =T (r) 

(2) 

T * - °r “e * “e 2 ) 2 

(3a) 

1 


r - ^ («r 2 +«r‘. +£ e 2 ) 2 

(3b) 

£ r = | r T(« r -|“e> 

(4a) 

£ « * 1 7 - i °r> 

(4b) 

S = t £ | K ♦ “*)] 

(4c) 


The constants l/2 and l/3 in equations (4) are determined from the 
condition defined by equation (l) and one of the equations (3). 
Only five of these equations are therefore independent. 
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The finite-strain concept (references 8 to 10), which con- 
siders the instantaneous dimensions of the element, is used because 
large deformation in the strain-hardening range is considered. 

The stress is then equal to the force divided by the instantaneous 
area and the strain to the logarithm of the instantaneous length 
divided by the initial length of the element (references 8 to 10 ) . 
It is mentioned in reference 5 that as long as the deformation the- 
ory is applicable, the logarithmic strain can also be used. The 
strain-displacement relations for this problem are then as follows: 


10g e 

d(r+u) 

dr 

(5a) 

l°g e 

r+u 

r 

(5b) 

3og e 

h 

(5c) 

^init 


From the condition of equilibrium in the radial direction of 
the small element in figure 1(b), the following equation of equi- 
librium is obtained (reference 5): 


(r+u) 


d(a r h) 

d(r+u) 


(°e - 0 r) h 


l 2 h r+u dr 
’ init r d(r+u) 


( 6 ) 


Nine equations defining this problem are equations (l), (2), (3b), 
(4a), (4b), (5a), (5b), (5c), and (6), which involve nine 
unknowns: cx r , Oq, £ r , € q, e z , y, 1 , h, and u. Equations (l) and 

(5) can be used to eliminate € z , u, and h, resulting in a com- 
patibility equation. The equations defining this problem are then 
reduced to six equations with six unknowns: a r , Oq, £ r , £q, y, 

and t. Two of the four unknowns, o r , Oq, £ ^, and Cq, may be 

eliminated by using equations (4a) and (4b) or (3b). The quan- 
tity T is a known function of y , which is experimentally deter- 
mined from a simple tensile test. The problem is then reduced to 
one involving three unknowns. Obtaining the solution of the 
resulting equations is not, however, a simple matter. This diffi- 
culty can be avoided by using the following transformation (refer- 
ences 1 and 5 to 7) because the yielding surface of plane plastic 
stress based on the deformation theory is an ellipse: 
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0 0 + a r = 3 / s/i" T sin a 
0g - CT r = 1 cos a 


Rewriting the principal stresses as functions of T 

®r - ’ <Vs sin a - cos a) 

T sin a + cos a) 


and a gives 



where T , which is a function of X, varies with r and also with 
load. With equations (7) substituted into equations (4a) and (4b), 
the principal strains can be expressed as functions of x and a 
as follows : 


e = — (sin a - f \Jz cos a) 

2V2 

€ q = — f — (sin a + V5~ cos a) 
y 2^/2 


( 8 ) 


In these equations, the parameter a is closely related to the 
ratio of principal stresses, inasmuch as 


a r sin a - cos a 

sin a + cos a 


and varies almost linearly with a for the range of a encountered 
in the present problem, which is the same range as in the thin 
plate with a circular hole as shown in figure 2 of reference 6. 

By using the transformation, the equations for this problem are 
reduced to the following two nonlinear differential equations, 
which are to be solved with an experimentally determined func- 
tion t(t) : 


A 


D 



+ 


+ 



C 


F 
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or 


/r\ da 

w HD = 


CE - FB 
AE - DB 


>v 


dr FA - CD 
EA - BD 


l (D 


0) 


J 


where 




. i— . . /— . r cos a 

A = [^J 3 cos a + sin a) - (V3 sin a - cos a) — 

, rz . S (r & T T Sin a\ 1 

= (Vo sxn a - cos a) ^ - 

--if » i "*) 2 £ (I) 


C = 2 (cos a) e 

D = (^3 sin a - cos a) y 
E = - (V3 cos a + sin a) 

(-\/t r °° s “ 


(S r cos °) rr , , ,_>2 V? 


sm a 


>( 9 ») 


= 2^” 


1 - e 


J 


f’(r^a-) 

Using equations (9a) and expanding the terms of e into a 

series result in 

> 

CE - BF = - 2HL - 2^3 HJ^f g J £ ) + ^(t)(|) f 2 (a ' r) 

AF-CD = | 8 H 2 - 2^~HL [l - f x (a,r)] + ^(t) f 2 ( a ^)J ^ 

1 1 £L _./l i r_ r ) 

\t dr v 2 j r 0 7 y 


)(9b) 


AE - BD = - L 2 - J 2 
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where 


fl(a.,r) = 


4 


(cos a)y 


1 - e 


(cos a)r 

= i - (cos a) y + \ (cos 2 a) y 2 . . . 


f 2 (a,y) = e 


7 ? 


sm a 


= 1 + 


(?? 8in " + 2 Sin “j r2 + Ut? Sln °j 


r + • • 


H = cos a 

J 
K 1 


- ^ sin a - cos a 
(a>k ) 2 




L cos a + sin 


a 


LINEARIZATION OF EQUATIONS 


In equations (9b), the variable t occurs not only in the 


combination of Jr t~ as in the case of a thin plate with a circu- 
T dy 

lar hole (reference 6), but also in the loading term. As proposed 

r dT 

in reference 6, the term j can be replaced by m, which is 


equal to the slope of a straight line approximating the T(y) curve 
on the logarithmic plot within the range of y encountered along 
the radius of the disk, y b to y Q . Thus m is a function of y Q j 

in other words, the value of m for one material is different for 
different loads. The special case where a straight line is used to 
approximate the whole strain-hardening range of the T(y) curve 
of a given material on the logarithmic plot, so that m is con- 
stant through the strain-hardening range, is the well-known power- 
law approximation. As in the case of a thin plate with a circular 
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hole, "the power-law approximation gives very good approximate 
results for this problem, as will be shown in the section "Calcu- 
lations, Results, and Discussion." The term - — can thus be 

t dr 

replaced by m without appreciable error. 

Furthermore, the term e< l ua ^ : ^ ons (9b) can be 

replaced by a constant C^, which is determined in a manner simi- 
lar to that in the case of a thin plate with a circular hole (ref- 
erence 6), as will be described. The terms (— — - A I— 2: X. y ) 

V dr V2 JL 7 

in equations (9b) then become (m - C-jr 0 ) • 

The general information concerning the effect of the 
T(r) curve of the material on the solution of this problem can be 
obtained only if the quantity t/t q in the loading term can pos- 
sibly be expressed as a function of (m - C-jj 0 ) . For any part of 

the T (r) curve that does not deviate greatly from a straight 
line in the strain-hardening range within the value of y con- 
sidered, r D to Tv the following equation can be written with 
sufficient accuracy: 



T - T C - (r 0 - r)Q 

av 

- T o 1 - (l 

'»> 

where 

(^j is the average slope of t(r) 

from r o to r b 

and 

(t q > t d ) is a known point 

on the T(y) 

curve . The quan- 

tity 

To ( d T \ 

=— I- 3 — I does not equal 

0 ' a ' ' av 

0 

O 

1 

£ 

but has a certain 


relation to it. Results obtained in reference 5 show that the val- 
ues of a and t/t o are not very sensitive to the 

of the material. If only a linear relation between 

and (m - C-jy 0 ) is retained, equation ( 10 ) becomes 

■f = 1 " t - f) C 2 (m ' C l r o> 

O \ O' 


t(t) curve 



( 10 a) 
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An approximate value of* C-^ = 0.5> which is equal to the mean value 
of »/— i — along the radius, can be used in equation (10a) . The 

\2 j r 0 

values of T 0 > 'f\ ) /Xo > an< ^ m are known ^ or point (t 0 > ' 

By using the values of t/ t q = t^/t q taken from the true octahed- 
ral shear stress -strain curves of several materials and using equa- 
tion (10a) and X /r Q = X\,/x 0 , the constant C 2 is then found to 

be 


C 2 = 2 


A zero value of C 2 is used for the case of negative (m - C-jj 0 ), 
because there is no negative slope on the true stress-strain curve 
of any material. Thus: 

= 1 - 2 f 1 - (m - C lTo ) (10b) 


It may be noted that this relation is very approximate and is only 
used to determine the variations of a and x/Xq along the radius, 

which are not sensitive to the t(t) . curve of the material. For 
the stress and rotating speed, which are dependent on the 
T(y) curve, this approximate relation cannot be used. The 
stresses are determined from the values of a and x/x Q obtained 
by the partly linearized solution and the true octahedral shear 
stress-strain curve of the material obtained from the tensile test. 
The rotating speed is determined from the distributions of tangen- 
tial stress and strain. 


Substituting equation (10b) into equations (9b) and the result- 

(x dT Id 1 X \ 

ing equation into equations (9), replacing ( — - yTj j — x) 

by (m - C-jTo), and neglecting the small terms (as in reference 6), 
the following equations are obtained: 
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da 

K!) 


<3 

(!) 


2H 

L 


1 + (m-Cq r c 

>> - I 1 

| ^ 

DO 

[i - 2| 

l 1 - 3 ( - ci r °>] 

r 


t 1 + © 

2 

(m-( 

: 1 r 0 )j 

l© 




8H 
T 2 


i . ^ (; 

l 6H 2 

!) 2 f 


(m-Cq r Q )j 

ll 

\- 
' To 


[l + (£) (m-c x r 0 )J 

l(!) 


>(H) 


J 


Let 


n = m - C 


1 


( 12 ) 


In equation (12), once the value of x Q is chosen, m can be 
determined from the t(t) curve of the material; Cq is a con- 
stant and n is therefore a constant determined from the 
T(f) curve of the material. Substituting equation (12) into equa- 
tions (ll) gives 

•\ 



i • 2 i 1 - 6 n ]' 1 © 2 

d (!) L [i + (# n ](D 

4 (3 3/ f 1 - M 1 • 


d (!) l2 [ 

d ♦ ©3(1) 


( 13 ) 


J 


In equations (13), all the terms are functions of a, x/x 0 > and 
r/k except Kq, which is an arbitrary loading constant, and n, 
which is also a constant during the calculation. The values of a 
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and t/t d c an be obtained along r/k until a reaches the value 
that satisfies the outer boundary condition (a = ^-) . Because r 

at the outer boundary is equal to b, the constant k can be 
determined. (For details, see reference 5.) The variations of a 
and t/t q with r/b for different values of n can then be cal- 
culated. Compare the values of a and Y"/t d along r/b obtained 

from equations (13) for several values of n with the values of a 
and t/t q along r/b obtained in reference 5. (The values of n 

for these last cases are calculated from equation (12) with the 
approximate value C-^ = 0.5.) The variations of a and r/r 0 

with r/b having the same n values obtained from equations (13) 
and from reference 5 are quite close; therefore = 0.5 can be 

used and equation (12) becomes 

n = m - 0.5 r o (12a) 

The relations of a and r/r Q with r/b for several values of n 

are given in both tabular and curve form to facilitate solution of 
this problem for any material under any maximum strain. This 
method of determining the distributions of a and t/t o along r/b 

will be referred to as the "partly linearized solution." 


PRINCIPAL STRESSES AND STRAINS 


After the variations of a and y/v 0 with r/b are obtained 

by the partly linearized solution, the principal stresses and 
strains can be obtained from equations (7) and (8) together with 
the actual t and y relation of a given material. The equations 
show that e r /y, e 0 /r, cr p /r , and ct 0 /t are functions of a only; 


they can be calculated for different values of a and given in 
tabular form to facilitate solution of any given case. Equa- 
tions (7) and (8) can be written as 




sin a 


cos a) 


sin a + cos a) 


(7a) 




MCA IN 2367 


13 



(8a) 


The values of principal stresses and strains can then be obtained 
by simple multiplication. 

It is usually interesting to know the ratios of maximum and 
minimum tangential stress and strain, because the rotating speed 
of the disk is mainly determined by the tangential stress and 
strain along the radius . The maximum tangential stress and strain 
occur at the center of the disk; at r = 0, 


The minimum tangential stress and strain occur at the rim of the 
disk; at r = b, 


From equations (7), the ratio of maximum to minimum tangential 
stress is 


From equations (8) , the ratio of maximum to minimum tangential 
strain is 





T 

O 


( ff0 )b T b ( sin | + cos |) 


( £ e) 0 r c ( sin | ^ cos |) To 
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IDEALLY PLASTIC AND APPROXIMATE SOLUTION 
For the case of n = 0, equations (13) become 



(13a) 


(13b) 


Equations (13a) and (13b) can be integrated numerically} these equa- 
tions are much simpler than those for the cases of n / 0. For 
ideally plastic material with the infinitesimal strain concept, 

r is infinitesimal and m = 0} the constant n is then equal 
o 

to zero} this case is a special one of n = 0. For this case, 

Nadai (reference l) obtained a similar relation between a and 
r/b as equation (13a) . 

Equations (15a) and (13b) can also be used to determine 
approximately the variations of a and y/t 0 with r/b for a 

material with n / 0. An approximate method of solution is there- 
fore proposed. The procedure of the approximate solution is as 
follows:' Use the variations of a and x/x Q with r/b for 

n = 0 (or for ideally plastic material with infinitesimal strain) 
as the approximate solution of the variations of <x and y/y o 

with r/b. The principal strains e g and e r can be calculated 
from the approximate relations of a, x/x q , and r/b because the 

strains are functions of ct and y only. The stresses, which 
are functions of a and t, can be calculated from equations (7a) 
by using the tensile true T(r) curve of the material and the 
relations of co, x/x 0 , and r/b determined for n = 0. 


DETERMINATION OF ROTATING SPEED OF DISK 

p 

The rotating-speed function p(cob) can he determined hy con- 
sidering a circular sector as sho^rn in figure l(c) . The radial 
component of force acting on the sector due to Qq is equal to 
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b 

„ , . A 9 d(r+u) , 

2a a h sin i dr 

0 2 dr 

The centrifugal force on the same sector is equal to 


P h init r ( A0 ) “ 2 ( r+u ) ^ 

From the condition of equilibrium of the sector in the radial 
direction, the following equation is obtained: 





The value of p(tob)' 
stress and strain. 


is then determined only by the tangential 
Accurate values of p(aob)^ can be obtained by 
substituting into equation (15) the tangential stress and strain, 
Oq and Cq , calculated from the linearized solution. 


APPROXIMATE VALUE OF ROTATING SPEED OF DISK DETERMINED 

DIRECTLY FROM TENSILE STRESS -STRAIN CURVE OF MATERIAL 

If only the approximate value of the rotating-speed func- 
tion p(cbb)^ is required, equation (15) can be made even simpler. 
Consider first the denominator because it is a function of &, 

and r/b only. Expanding e into a series, the denominator 
becomes 
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where ir 0 is a constant and can he taken out of the integral sign. 
Thus the denominator of equation (15) becomes 


VJ 


pi 

£ e 

e 

0 



= 3 + °1 r c 


+ °2 ^o 



(16a) 


where G^, G£, and so forth are integrals, which sire functions of a, 
x/x q > and r/b only. The integrals are given in table I. 

The integrals can be calculated using the values of a and 
x/x o along r/b obtained from the approximate solution. The var- 
iations of a and t/t d with r/b are independent of material 
according to the approximate solution. 

The numerator of equation (15) is a function of and 

r/b. The tangential strain Cq can be treated in a manner similar 
to that used for the denominator. The tangential stress Oq can 
be written as where a^/ T is a function of a, and T 

can be written as an approximate function of y, thus 

where m is the slope of a straight line approximating the tensile 
T(y) curve of the material on the logarithmic plot within the 
range of r encountered along the radius of the disk Tq, to r 0 , 

and K is a constant that can be determined from the actual 
T (y ) curve within the same range of y . The integrand of the 
numerator thus becomes: 
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Expanding the exponential function of e into a series gives 

r Q 

-(m logg — + c Q ) 

e 

t/t ^o , NlI/t V »2 1 , ^o .>3 

= 1 - (m log e — +e e ) (m log e —+e 0 ) - jr (m log e — + +• • • 



The values of m and r 0 remain constant along the radius and can 
therefore he taken out of the integral sign; thus the numerator of 
equation (15) becomes 



" I 4^o + I 5 m T 0 + • • •) (16b) 

The integrals I Q , I-j_, I 2 , I 3 , 1^, Ig, and so forth are functions 
of a, y/y 0 > and r/b, and are independent of material according 

to the approximate solution. These integrals are also given in 
table I. Substituting equations (16a) and (16b) into equation (15) 
and neglecting the small terms result in 

p(cob) 2 _ Ip - Iqm + Igm 2 - I 5 m° - I 4 r 0 + I 5 mr 0 

| + G l r o + G 2f o 


T 

O 


(17) 
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The values of I Q , I , . . . and G 2 are calculated from the 

formulas listed in table I by using the a, r/r Q , and r/b rela- 
tions obtained in the case of n = 0 (the approximate solution) . 
These values are as follows: 


2.3031 

*4 

0.5239 

2.0281 

*5 

.3514 

1.3131 

G 1 

.0537 

.6362 

°2 

.0047 


Equation (17) shows that p(cbb) 2 / T 0 is a function of m and r o - 
The value of y Q is the maximum value of r in the disk. Once 
the value of r Q is selected, m can be obtained from the tensile 

r(y) curve of the material within the range of r considered 
(r^ to r Q ; Tb/r 0 is approximately equal to 0.15), and is 

the value of t corresponding to r o taken from the T (y) curve 
of the material. Equation (17) thus gives the approximate value 

of p(cob) 2 , which is directly determined from the tensile stress - 
strain curve of the material. Curves of p(a5b) 2 /T Q against y Q 

for different constant values of m can be calculated and plotted. 
The value of p((bb) 2 /T 0 can then be read directly from the chart 
for given values of m and ir 0 . 

CORRELATION BETWEEN EXPERIMENTAL BURSTING SPEED AND 

CALCULATED VALUE DETERMINED DIRECTLY FROM 

TENSILE STRESS -STRAIN CURVE OF MATERIAL 

It is interesting to compare the experimentally determined 
rotating speed corresponding to a value of at the center of 

the disk with the rotating speed predicted from the chart (or equa- 
tion (17)) mentioned in the previous section. In ordbr to deter- 
mine these speeds experimentally, it would be necessary to measure 

r at the centers of the rotating disks during operation. Such 
o 

data are unavailable in the literature to the best knowledge of 
the author. If it is assumed that fracture in the rotating disk 
takes place at the same value of r Q as the fracture octahedral 

shear strain in simple tension, then the data in reference 4, 
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together with additional data obtained at the Lewis laboratory for 
the work reported in reference 4, are sufficient to allow comparison 
between prediction and experiment . 

The predicted rotating speed is determined directly from the 
tensile octahedral shear stress -strain curve of the same material 
as the disk. These curves are calculated from the unpublished 
tensile test data for reference 4. The tensile specimens were cut 
in the radial direction of the disk of the same material and same 
heat treatment. The specimen was 3.00 inches long with a 1.00-inch 
gage length. Detailed dimensions of the tensile specimen are 
given in figure 4 of reference 4. During the tensile test, the 
load-elongation curve was obtained by a recording extensiometer. 

The original and breaking diameters of specimens were measured 
before and after the experiment. The load at which the. specimen 
broke was read from the dial of the tensile machine. From these 
data, a T(y) curve of the material can be obtained. First, cal- 
culate t and y by using the values taken from the load- 
elongation curve up to the maximum load point: 


r = 'Jz &L =*Jz log e (l+AZ) 

t = sA F _ p = Vf = _l 

3 a i 3 A 3 A 0 A SAq 


(1+ Al) 


(18) 


where 

A l elongation of 1-inch gage length 

A instantaneous area 

Aq original area 
l instantaneous length 


l Q original length 

P load of tensile specimen 

The t and x relation cannot be obtained from the tensile 
load-elongation curve after the maximum load is reached because 
of the nonuniform elongation within the gage length. However, the 
t and x fracture can be calculated from the original cross- 
sectional area, fracture area, and fracture load: 
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T fr 


P fr 
3 A fr 


fr 


r- A- 

" lo ®e — 
rr 




> 


(18a) 


where the subscript fr means at fracture. Plot the t and y 
relation obtained on a logarithmic scale. Draw a straight line 
to approximate the T(y) curve in the range of ) to T fr ; 

for a solid rotating disk, y^y is approximately equal to 0.15 
(for the case of n = 0) . The slope of this straight line gives 
the value of m. Find the value of p(oob) /t corresponding to m 
and from the chart. The predicted value of p(tob) when 

the disk bursts will be equal to £p(<x>b) /t 1 T fr by assuming the 
disk breaks when the maximum y in the disk (y at the center 
of the disk) reaches the value of y^ r in simple tension. 


CALCULATIONS, RESULTS, AND DISCUSSION 

In order to observe the degree of the approximation resulting 
from the use of a power function representing the octahedral shear 
stress-strain relation over the complete strain-hardening range 
(power -law approximation) , a calculation is made with the follow- 
ing relation between t and y: 

t ( r) = 126,000 r 0 ' 25 

The - constants are chosen to approximate the T(y) curve of 
Inconel X over the whole strain-hardening range, as shown in fig- 
ure 2. In figure 3, the variations of a, y> cr r , Oq, e r , and Cq 

with r/h for the case of T 0 = 0-30 obtained by the power-law 

approximation are compared with the values obtained from the actual 
r(y) curve of Inconel X; it can be seen that the power-law 

approximation gives good results. These results indicate that 
replacing the term V equations (9b) by m should introduce 

very little error. Also, if a simple analytical function of t (y) 
is desired for analysis, the simple power function representing 
the t and y relation of a given material will give a very good 
approximation . 
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The values of a and x/x Q along the radius are calculated 

for n = -0.05, 0, 0.1, 0.2, and 0.3 from the partly linearized 
equations (13) . The variations of a and x/x Q with r/b for 

the different values of n are plotted in figure 4; these values 
are given in table II. The values of a and x/x 0 along r/b 
obtained in reference 5 for 16-25-6 and Inconel X are plotted 
in figure 5; the values of n for these cases are calculated from 
equation (12a) and indicated on these curves. The curves obtained 
from the partly linearized solution in figure 4 having the same 
range of values of n as those obtained from reference 5 are also 
plotted in figure 5 for comparison. The curves for the same 
n value of the two solutions agree very well for most cases. For 
simplicity, the solution based on the deformation theory obtained 
in reference 5 is designated the exact solution. The simple rela- 
tion n = m - 0.5 x 0 can then be used as a good approximate cri- 
terion to find the variation of a and x/x Q with r/b for dif- 
ferent materials and different maximum strains. This generality 
leads to a general criterion of applicability of the deformation 
theory of plasticity to this problem for any material in the strain- 
hardening range. The criterion is that if the value of n for a 
given material is constant or approximately constant in the strain- 
hardening range, the deformation theory can be applied to this 
problem for the material. For the special case of infinitesimal 
strain, the condition of n being constant reduces to m being 
constant, which is the same condition obtained by Ilyushin 
(reference 11) . 

In figure 4(a), the maximum variation of a is about 15 per- 
cent whereas n varies from -0.05 to 0.3. For most materials, m 
increases with X } thus it can be seen from equation (12a) that 
m and x affect n in an opposite sense, so that n does not 
change much with strain. The value of n for most materials of 
any maximum strain r 0 varies from -0.05 to 0.3 (references 12 

to 14); only a small part of this variation is due to the varia- 
tion with strain for a given material. The variation of a during 
loading is expected to be very small for most materials; conse- 
quently the deformation theory of plasticity is applicable to this 
problem within engineering accuracy. 

The values of cr r / t, o Q fr , C r /r, and c q/x are calculated 

from equations (7a) and (8a) for the range of a considered and 
are given in table III. The distributions of principal stresses 
and strains are obtained by the partly linearized solution for the 
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cases of r Q = 0.3000 for Inconel X and r Q = 0.4400 for 16-25-6. 

These results are plotted in figure 6. The results obtained by 

the exact solution are also plotted for comparison. The examples 

with r D = 0.3000 for Inconel X and r Q = 0.4400 for 16-25-6, with ^ 

the values of n equal to 0.15 and -0.05, respectively, are V 

chosen because they give the largest deviation of the examples con- 
sidered between the partly linearized and exact solutions in curves 
of a and t/t q against r/b. For y Q = 0.3000, Inconel X, 

n is equal to 0.15 and the values of a and x/x Q corresponding 

to different values of r/b can be obtained from table I by lin- 
ear interpolation. In figure 6, the principal stress and strain 
distributions obtained by the partly linearized solution are shown 
to agree very well with those obtained by the exact solution, even 
for these two cases that have the largest deviation of the varia- 
tions of a and t/t c with r/b between the two solutions. For 

the rest of the cases, the a and r/r 0 curves obtained by the 

two solutions are approximately the samej therefore the distribu- 
tions of principal stresses and strains will also be practically 
the same. 

The variations of a and t/t 0 with r/b obtained by the 

approximate solution are the same as those for the case of n = 0 
in figure 4 . The principal strains obtained by the approximate 
solution and the exact solution are plotted in figures 7(a) and 
7(b) for 16-25-6 and Inconel X, respectively. The stress distri- 
butions obtained by these two solutions are plotted in figures 7(c) 
and 7(d) for 16-25-6 and Inconel X, respectively. The stresses 
and the strains obtained by the approximate solution give fairly 
good results . 

The variations of a, x/j Q ) q/ {*■ q) 0 > and Sr/^r^o with 

radius for ideally plastic material are the same as those obtained 
by the approximate method. The principal stresses for the ideally 
plastic material with y Q equal to 0.4400 and 0.3000 are also cal- 
culated. In this case, T is constant and is chosen equal to the 

t for which the solutions are being compared. These distribu- 
o 

tions of principal stresses along the radius as well as those 
obtained by the partly linearized, exact, and approximate solutions 
are plotted in figures 8(a) and 8(b) for the case of r 0 = 0.4400 

for 16-25-6, and r Q = 0.3000 for Inconel X, respectively, for com- 
parison. It is seen from these figures that the stresses (espec- 
ially the tangential stress, which is the most important part of 
the solution) obtained by the ideally plastic material are not good 
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enough to represent the actual results. The curves obtained by 
the three other methods agree very well. These results also indi- 
cate, as in the case of the thin plate with a circular hole (ref- 
erence 6) that the ratios of strains along the radius to the maxi- 
mum value and the ratio of the principal stresses are essentially 
independent of the t(t) curve of the material and the maximum 
strain of the disk, but that the stresses are very much dependent 
on the material. The conclusions obtained in reference 5 for the 
problem of the rotating disk (with body force) of Inconel X and 
16-25-6 are then expended to most materials. 

2 

The rotating-speed function p(o)b) is calculated from equa- 
tion (15) by using the values of Og and £g obtained in the 

partly linearized solution for Inconel X and 16-25-6. These values 

of p(cdb) and those obtained from the exact solution are plotted 
against r 0 in figure 9. The loads obtained by these two methods 

are practically the same (the two curves coincide) . The relations 
between the approximate value of p (cub ) 2 / T q and r Q for several 

values of m are calculated from equation (17) and plotted in fig- 
ure 10. The variations of load p(o>b) 2 with r Q for Inconel X 

and 16-25-6 are obtained from figure 10 and are plotted in fig- 
ure 9 for comparison. The values obtained by the approximate sol- 
ution are also very close to the values obtained by the exact 
solution. The percentage deviation of p(<nb) 2 of the approximate 
solution from the exact solution is much less than the percentage 
deviation of tangential stress and strain. This can be explained 

p 

by the fact that p(oib) (equation (15)) increases with increas- 
ing value of ag but decreases with increasing value of €g} 
values of Og and e g obtained by the approximate solution 

(fig. 7) differ from those of the exact solution in the same 
direction; hence the part of error cancels. The percentage devia- 
tions of radial stress and strain in some cases are quite large, 
but fortunately these quantities do not enter in the calculation 
of p (db ) u . 

2 

Use is made of figure 10 to calculate p(cob) from the tensile 
stress -strain curve of the same materials as the solid disks con- 
sidered in reference 4. (in this reference, data were obtained 
for disks that showed no imperfections upon X-ray and surface 
inspection and for disks in which defects did exist. For the 
present purposes, only the disks having no defects were considered.) 
From the load-elongation curve, original diameter, final diameter, 
and fracture load in tension, which data were obtained at this lab- 
oratory for the work reported in reference 4, the t (y ) curves of 
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these materials were calculated using equations (18) and (18a) and 
are plotted in figure 11 on a logarithmic scale. Six materials 
(SAE 1078 steel of two different heat treatments, beryllium- copper 
alloy of three different heat treatments, and one nickel-base 
alloy) are considered. For simplicity, these materials are herein 
designated materials A, B, C, D, E, and F. The values of T|. r 

and Tf x f° r these materials are given in table IV; the reduction 

in nonnecked diameter, which was used as a parameter in some fig- 
ures of reference 4, is also given in table IV in order to corre- 
late the results given in reference 4 and herein. (This reduction 
in nonnecked diameter is the reduction in diameter in the region 
of uniform strain of the specimen.) In table IV are listed the 

p 

values of m and p(oib) , which are calculated by taking r Q at 
the center of the disk equal to Tf r of the tensile test, except 

for material C. Because the load-elongation curves of the tensile 
test show that material C has a large necking effect, several val- 
ues of p(o)b) 2 for different y Q between y = 0.1342 (the y 
at the maximum load of the tensile test) and y fr were calculated 

and are given in table V. It can be seen that the calculated 
rotating-speed function p(cub) 2 increases to a maxim-urn at 
y^ = 0.3 to 0.4 and then decreases. This fact indicates a case 

of instability of the rotating disk. 


The experimentally determined values of p(cbb) 2 are also 
given in table IV. The percentage difference of calculated p(cbb) 2 
to the experimentally determined p (cob ) 2 as well as the percentage 
difference of calculated cjd to the experimental ly determined o 
is given in the same table. It can be seen from these values that 
the percentage difference in cn .is small. More reliable calcu- 
lated values of p(oab) 2 may be obtained, however, if: (l) the 

T(y) curves obtained from the tensile test are corrected for the 
triaxiality and nonuniform stress distribution introduced by neck- 
ing, and (2) the T(y) curve of a few tangential tensile specimens 
taken at different radii of the disk (including a specimen passing 
through the center of the disk) are used, because the material of 
the disk is not quite uniform and p((»b) is a function of £q , 

Oq, and r/b, as shown in equation (15) . The calculated values of 

p 

p (oib) at disk failure are plotted against the experimentally deter- 

mined p(aob) 2 in figure 12. It can be seen that these points are 
slightly above the 45° straight line. Calculations based on 
t(t) curves of tangential tensile specimens may give lower values. 
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It should be emphasized that a! 1 the results and discussion 
are only true for this problem tinder plastic deformation in the 
strain-hardening range in which the elastic strains are small com- 
pared with the plastic strains and there are neither time and tem- 
perature effects nor unloading. 


CONCLUSIONS 

The results obtained for the rotating disk in the strain- 
hardening range, in which the elastic strains are small compared 
with the plastic strains, lead to the following conclusions, which 
are similar to those obtained in the case of the thin plate with a 
circular hole: 

1. The results obtained by the partly linearized solution 
agree very well with those obtained by the exact solution based on 
the deformation theory of plasticity. The amount of computation 
is greatly reduced by the partial linearization. The stresses and 
the strains can be obtained for any material under any maximum 
strain by a simple multiplication using the tables or curves given 
in this paper. 

2. The variation of a parameter, which is determined from the 
tensile-stress-and-strain curve of the material, can be used as a 
simple general criterion of the applicability of deformation the- 
ory for this problem. 

3. The results previously obtained for Inconel X and 16-25-6, 
namely, the variation of the ratio of the principal stresses with 
radius and the ratios of the strain along the radius to their maxi- 
mum value are essentially independent of the octahedral shear 
stress -strain relation of the material but that the distributions 
of the stresses depend very much on the material, are extended 

to most materials. 

4., Results obtained for the ideally plastic material with the 
infinitesimal strain concept give good approximate values of prin- 
cipal strains but not principal stresses. 

5. Sufficiently accurate values of principal stresses can be 
obtained by an approximate method in which the stresses are cal- 
culated by using the strains obtained from the ideally plastic 
material together with the actual octahedral shear stress -strain 
relation of the material. 
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6. The rotating speeds of the disk obtained by the partly lin- 
earized solution are practically the same as those obtained by the 
exact solution; and those obtained by the approximate solution, in 
which the rotating speeds are determined directly from the tensile 
stress -strain curve of the material, are also very close to the 
exact value. 

7 . Good correlation between experimental bursting speed and 
the calculated value determined directly from the tensile stress- 
strain curve of the material is obtained. 

8. If a simple analytical function representing the octahedral 
shear stress-strain relation is required for analysis, the power- 
law approximation can be used. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 23, 1951. 
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TABLE I 


- VALUES OF I ID G USED IN EQUATIONS (16) AND (17) 
PI 
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TABLE II - VARIATION OF a AND r/r o WITH r/b 
FOR VARIOUS VALUES OF n 


a 

n = - 

0.05 

n = 

0 

n = 

0.1 

n = 

0.2 

n = 

0.3 


r/b 

r/r 0 

r/b 

r/r D 

r/b 

r/r D 

r/b 

r/r 0 

r/b 

r/r c 

1.5708 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

1.5675 

.0480 

.9850 

.0510 

.9850 

.0600 

.9850 

.0690 

.9850 

.0800 

.9850 

1.5615 

• 0750 

.9675 

.0900 

.9675 

.1010 

.9675 

.1140 

.9675 

.1260 

.9675 

1.5480 

.1170 

.9250 

.1325 

.9250 

.1565 

.9250 

.1780 

.9250 

.2050 

.9250 

1.5300 

.1555 

.8700 

.1780 

.8710 

.2050 

.8730 

.2350 

.8750 

.2720 

.8770 

1.5045 

.2025 

.8010 

.2250 

.8030 

.2580 

.8060 

.2970 

.8090 

.3380 

.8120 

1.4820 

.2370 

.7490 

.2600 

.7500 

.2970 

.7550 

.3400 

.7590 

.3840 

.7630 

1.4490 

.2820 

.6780 

.3050 

.6820 

.3475 

.6880 

.3930 

.6920 

.4400 

.6960 

1.4105 

.3280 

.6100 

.3510 

.6150 

.5975 

.6210 

.4420 

.6280 

.4910 

.6350 

1.3670 

.3775 

.5420 

.4000 

.5490 

.4460 

.5580 

.4920 

.5675 

.5410 

.5760 

1.3180 

.4280 

.4790 

.4505 

.4875 

.4950 

.4995 

.5400 

.5120 

.5900 

.5210 

1.2635 

.4820 

.4200 

.5025 

.4295 

.5450 

.4450 

.5890 

.4575 

.6370 

.4710 

1.2045 

.5350 

.3690 

.5540 

.3790 

.5940 

.3950 

.6360 

.4120 

.6800 

.4250 

1.1400 

.5880 

.3230 

.6060 

.3345 

.6440 

.3520 

.6830 

.3695 

.7230 

.3850 

1 .0705 

.6440 

.2840 

.6600 

.2950 

.6930 

.3145 

.7290 

.3325 

.7650 

.5500 

.9965 

.7000 

.2505 

.7150 

.2610 

.7440 

.2820 

.7730 

.3000 

.8050 

.3175 

.9175 

.7575 

.2220 

.7700 

.2330 

.7930 

.2525 

.8170 

.2720 

.8440 

.2915 

.8330 

.8140 

.1975 

.8240 

.2075 

.8420 

.2280 

.8610 

.2480 

.8820 

.2690 

.7430 

.8710 

.1760 

.8780 

.1870 

.8920 

.2070 

.9040 

.2270 

.9180 

.2490 

.6465 

.9290 

.1590 

.9330 

.1690 

.9410 

.1900 

.9490 

.2100 

.9560 

.2310 

.5236 

1.0000 

.1430 

1.0000 

.1520 

1.0000 

.1720 

1.0000 

.1930 

1.0000 

.2150 
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TABLE III - VALUES OF , o Q f T, cjr , AND £ 0 /r 

FOR CORRESPONDING VALUE OF a 


a 

a r/ T 

V T 

£ r/r 

t-e/r 

1.5703 

2.1213 

2.1213 

0.3536 

0.3536 

1.5675 

2.1172 

2.1253 

.3515 

.3556 

1.5615 

2.1098 

2.1326 

.3478 

.3592 

1.5430 

2.0929 

2.1487 

.3395 

.3674 

1.5300 

2 .0696 

2.1695 

.3283 

.3782 

1.5045 

2.0355 

2.1978 

.3122 

.3933 

1.4820 

2.0044 

2.2216 

.2979 

.4065 

1.4490 

1.9568 

2.2544 

.2765 

.4253 

1.4105 

1.8936 

2.2896 

.2513 

.4468 

1.5670 

1.8295 

2.3253 

.2223 

.4702 

1.3180 

1.7476 

2.36G2 

.1892 

.4955 

1 . 2655 

1.6515 

2.3924 

.1518 

.5222 

1.2045 

1.5419 

2.4192 

.1108 

.5494 

1.1400 

1.4160 

2.4389 

.0655 

.5770 

1.0705 

1.2738 

2.4488 

.0165 

.6040 

.9965 

1.1157 

2.4463 

- .0358 

.6295 

.9175 

.9401 

2.4289 

- .0915 

.6530 

.8330 

.7458 

2.3935 

-.1503 

.6735 

.7430 

.5331 

2.3370 

-.2118 

.6902 

.6465 

.3003 

2 . 2555 

-.2758 

.7018 

.5236 

.0000 

2.1213 

- .3536 

.7071 


TABLE IV - DISK MATERIALS, TENSILE -TEST PROPERTIES, CALCULATED AND EXPERIMENTALLY 
DETERMINED BURST SPEEDS OF ROTATING DISK 







Burst speed 




Tensile-test 

properties 

Calculated from 
tensile t (y) curve 

Experimentally 

determined 

Difference 
of calcu- 

Difference 
of calcu- 

Material 8, 

Reduction in 
nonnecked 
diameter 
(percent) 

T fr 

T fr 

m 

(in range 
of 

r fr< r tAo> 

to r fr ) 

p(aib) 2 /r 0 

from 
fig. 10 

Calculated 
, .2 
p(a)b) 

Average 

calculated 

p(cob) 2 

p(u*) 2 

lated to 
experimental 

p(cub) 2 

(percent) 

lated to 
experimental 

CO 

(percent) 

Steel, 
SAE 1078 











A 

2.0 

68,900 

65,800 

0.1237 

.0897 

0.186 

.170 

5.61 

5.77 

386 X 10 5 
379 

382.5 X 10 3 

381.5 x 10 3 

0.26 

0.13 

B 

2.5 

52,500 

58,000 

.1298 

.2133 

.155 

.175 

5.75 

5.47 

301 

317 

309 

289.1 

292.8 

6.54 

5.19 

3.22 

2.56 

Beryllium 

copper 

C 

6.0 

78,120 

.8224 

.192 

4.15 

324 


310.2 

4.51 

2.16 

D 

9.0 

76,300 

1.3580 

.413 

2.72 

212 


206.2 

2.81 

1.40 

E 

16.9 

76,000 

1.3570 

.480 

2.60 

197.5 


191.4 

3.19 

1.57 

Nickel- 

base 

alloy 

F 

8.2 

94,700 

89,500 

0.2957 

.2538 

.305 

.294 

4.76 

4.92 

450 

440 

445 

433.9 

428.2 

2.56 

3.92 

1.34 

2.02 


a Two different heat treatments for SAE 1078 and three different heat treatments for "beryllium- copper alloy. 
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TABLE V - RELATION BETWEEN ROTATING SPEED p(oob) 2 
AND MAXIMUM OCTAHEDRAL SHEAR STRAIN AT CENTER 
OF DISK FOR MATERIAL C (BERYLLIUM-COPPER) 



T o 

m 

p(a>b) 2 /r 0 
from fig. 10 

Calculated 

p(oib) 2 

0.1342 

54,560 

0.092 

6.04 

330 x 10 3 

.2000 

58,000 

.112 

5.80 

336 

.3000 

62,800 

.152 

5.38 

337 

.4000 

66,200 

.162 

5.10 

337 

.5000 

69,800 

.185 

4.81 

336 

.8224 

78,120 

.192 

4.15 

324 
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Figure 1. - Rotating disk, an element, and a sector in deformed state. 
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Natural octahedral shear strain, 7 
(a) Linear scale plot. 
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Figure 2. - Octahedral shear stress-strain curves. 



(b) Logrithmic scale plot. 
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Figure 3® - Comparison of results obtained for t(y) for Inconel X 
and power-law approximation, t(y) = 126,000 yO-25 for rotating 
disk. y q = 0.30. 
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(b) Variation of octahedral shear strain with proportionate 

radial distance. 


Figure 3. - Continued. Comparison of results obtained for t(y) 

for Inconel X and power-law approximation, -r(y) = 126,000 y 0 * 25 
for rotating disk. Yo = 0*30. 
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Proportionate radial distance, r/b 

(c) Variation of principal stresses with proportionate radial distance. 


Figure 3. - Continued. Comparison of results obtained for t(y) for 

Inconel X and power-law approximation, t (y) = 126,000 ^0.25 for 
rotating disk. y 0 = 0.30. 
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(d) Variation of principal strains with proportionate radial distance. 


Figure 3. - Concluded. Comparison of results obtained for r{y) for 
Inconel X and power-law approximation, t(y) = 126,000 y°' 2 ~ for 
rotating disk. y q = 0.30. 
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Figure 4. - Variations of parameter a and ratio of octahedral shear strain to octahedral shear strain at center 
y/y Q with proportionate radial distance r/b obtained by partly linearized solution for different values of n. 
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Figure 4. - Concluded. Variations of parameter a 
center y/y 0 with proportionate radial distance 


and ratio of octahedral shear strain to octahedral shear strain at 
r/b obtained by partly linearized solution for different values of n. 
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Figure 5. - Comparison of variations of parameter a and ratio of octahedral shear strain to octahedral shear strain 
at center 7 /y 0 with proportionate radial distance r/b obtained by partly linearized solution and by exact solution 
for 16-25-6 and Inconel X. 
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Figure 5. - Continued. Comparison of variations of parameter a and ratio of octahedral shear strain to octahedral 
shear strain at center y/y Q with proportionate radial distance r/b obtained by partly linearized solution and by 

exact solution for 16-25-6 and Inconel X. 
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(c) Variation of ratio of octahedral shear strain to octahedral shear strain at center y/y Q with proportionate radial 

distance for 16-25-6. 


Figure 5. - Continued. Comparison of variations of parameter a and- ratio of octahedral shear strain to octahedral 
shear strain at center y/y Q with proportionate radial distance r/b obtained by partly linearized solution and by 
exact solution for 16-25-6 and Inconel X. 
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Figure 5. - Concluded. Comparison of variations of parameter a and ratio of octahedral shear strain to octahedral 
shear strain at center y/y Q with proportionate radial distance r/b obtained by partly linearized solution and by 
exact solution for 16-25-6 and Inconel X. 


. 


S2T2 


MCA TN 2367 


45 


I 1 

CO 



Figure 6. - Comparison of principal strains and stresses obtained by partly linearized 
solution and exact solution. 
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Figure 7. - Comparison of principal strains and stresses obtained by approximate solution and exact solution. 
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Principal stresses, o and o r , lb/sq In, 
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Figure 8. - Comparison of principal stresses obtained by different methods. 
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(b) Distributions of principal stresses along proportionate radial distance for y 0 - 0.5000, Inconel X. 
Figure 8. - Concluded. Comparison of principal stresses obtained by different methods. 
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Figure 9. - Comparison of relations between rotating -speed function 
p ( tub ) 2 and maximum octahedral shear strain at center of disk 
obtained by exact, partly linearized, and approximate solutions 
for Inconel X and 16-25-6. 
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Figure 10. - Variation of ratio of rotating-speed function to maximum octahedral shear stress at center of disk 
p(<i*) 2 /t 0 with maximum octahedral shear strain at center of disk for different values of parameter m. 

(A 14- by 15- in. print of this fig. is attached.) 
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Natural octahedral shear strain, 7 

Figure 11. - Octahedral shear stress-strain curves of tensile specimens cut from disks. 
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Calculated bursting-speed function, p(cub)2 
lb/sq in. 
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Experimental bursting- speed function, p(cnb)-', lb/sq in. 

Figure 12. - Relation between calculated and experimental bursting 

speed function. 
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